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ABSTRACT 

We report an analysis of the interstellar 7-ray emission in the third Galactic quadrant mea- 
sured by the Fermi Large Area Telescope. The window encompassing the Galactic plane from 
longitude 210° to 250° has kinematically well-defined segments of the Local and the Perseus arms, 
suitable to study the cosmic-ray densities across the outer Galaxy. We measure no large gradient 
with Galactocentric distance of the 7-ray emissivities per interstellar H atom over the regions 
sampled in this study. The gradient depends, however, on the optical depth correction applied 
to derive the H 1 column densities. No significant variations are found in the interstellar spectra 
in the outer Galaxy, indicating similar shapes of the cosmic-ray spectrum up to the Perseus arm 
for particles with GeV to tens of GeV energies. The emissivity as a function of Galactocentric 
radius does not show a large enhancement in the spiral arms with respect to the interarm region. 
The measured emissivity gradient is flatter than expectations based on a cosmic-ray propagation 
model using the radial distribution of supernova remnants and uniform diffusion properties. In 
this context, observations require a larger halo size and/or a flatter CR source distribution than 
usually assumed. The molecular mass calibrating ratio, Xco = N(H.2)/Wco, is found to be 
(2.08 ± 0.11) x 10 20 enr 2 (K km s" 1 )" 1 in the Local-arm clouds and is not significantly sensi- 
tive to the choice of H 1 spin temperature. No significant variations are found for clouds in the 
interarm region. 

Subject headings: cosmic rays - gamma rays: ISM - ISM: general 
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1. Introduction 

Knowledge of the distribution of cosmic-ray 
(CR) densities within our Galaxy is a key to un- 
derstanding their origin and propagation. High- 
energy CRs interact with the gas in the interstellar 
medium (ISM) or the interstellar radiation field, 
and produce 7-rays via nucleon-nucleon interac- 
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tions, electron Brcmsstrahlung and inverse Comp- 
ton (IC) scattering. Since the ISM is transpar- 
ent to these 7-rays, we can probe CRs in the lo- 
cal ISM, beyond direct measurements performed 
in the solar system, as well as in remote loca- 
tions of the Galaxy. Although mu ch effort has 



been made since the COS-B er a (e.g.. lStrong et al 



19881 : IStrong fc Mattoxl Il996l ). the results have 
been limited by the angular resolution, effective 
area and energy coverage of the instruments. The 
advent of the Fermi Gamma-ray Space Telescope 
enables studying the spectral and spatial distribu- 
tion of diffuse 7-rays and CRs with unprecedented 
sensitivity. 

Here we report an analysis of diffuse 7-ray emis- 
sion observed in the third Galactic quadrant. The 
window with Galactic longitude 210° < I < 250° 
and latitude —15° < b < +20° hosts kincmati- 
cally well-defined segments of the Local and the 
Perseus spiral arms and is one of the best regions 
to study the CR density distribution across the 
outer Galaxy. The regio n has been already stud- 
ied bv lDigel et all (|200ll ) using EGRET data. The 
improved sensitivity and angular re solution of the 



Ferm i LAT (Large Area Telescope; lAtwood et al 



20091) and recent developments in the study of 
the ISM allow us to examine the CR spectra and 
density distribution with better accuracy. We 
exclude from the analysis the region of the Mono- 
ceros R2 giant molecular cloud and the Southern 
Filament of the O rion-Monoceros complex (e.g., 



Wilson et all 120051) . in / < 222° and b < -6°, 



because 1) star forming activity and possible high 
magnetic fie lds suggested by the filamentary struc- 
ture (e.g., 



Morris, Montani. fc Thaddeus 
198' 



19801: 



Maddalena et al.l 119860 could indicate a special 
CR environment, and 2) an OB association in 
Monoceros R2 may hamper the determination of 
ISM densities from dust tracers (see § 2.1.2 for 
details). 

Study of the Aco conversion factor which 
transforms the integrated intensity of the 2.6 mm 
line of carbon monoxide, W^cOi into the molecular 
hydrogen column density, A(H2), is also possible 
since the region contains well-known molecular 
complexes. In the Local arm we find the molec- 
ular clouds associated with Canis Major OB 1, 
NGC 2348 and NGC 2632 (iMel'nik fc Efremov 
19951 : iKaltcheva fc Hilditchl l2000h . At a few 
kpc from the Solar System, in the interarm, 



lower-density region located between the Local 
and Perseus arms, we find M addalena's cloud 
( Maddalena fc Thaddeusl[l985h . a giant molecu- 
lar cloud remarkable for its lack of star formation, 
and the cloud associated wi th Canis Major OB 2 
(|Kaltcheva fc Hilditchl [2000h . 

This study complements the Fermi LAT study 
of the Cassiopeia and C epheus region in the second 
quadrant reported by lAbdo etaL ( 2010al ). The 
paper is organized as follows. We describe the 
model preparation in § 2 and the 7-ray observa- 
tions, data selection and the analysis procedure 
in § 3. The results are presented in § 4, where we 
also discuss the emissivity profile measured for the 
atomic gas and we compare it with predictions by 
a CR propagation model. A summary of the study 
is given in § 5. 

2. Modeling the Gamma-Ray Emission 
2.1. Interstellar Gas 

2.1.1. Hi and CO 

In order to derive the 7-ray emissivities as- 
sociated with the different components of the 
ISM we need to determine the interstellar gas 
column densities separately for each region and 
gas phase. For atomic hydrogen we used the 
Leiden / Argent i ne /Bo nn Galactic Hi survey by 



Kalberla et al. ( 20051 ). In order to turn the Hi 
line intensities into A(Hi) column densities, a 
uniform spin temperature Tg = 125 K has of- 
ten been adopted in previous studies. We will 
consider this option to directly compare our re- 
sults with th e former EGRET analysis of the 
same region (jDigel et al. 1 120011) and other stud- 
i es of the Galactic diffuse emission by the LAT 
(|Abdo et al.l l2009ai l2010ah Re cent Hi absorp- 
tion studies ( Dickev et al.l 2009h . however, point 
to larger average spin temperatures in the outer 
Galaxy, so we have tried different choices of Tg to 
evaluate how the optical depth correction affects 
the results. We will find that the emissivity per 
H 1 atom and the inferred CR density is affected 
by up to ~ 50% in the Perseus arm, and will take 
this uncertainty into account in the discussion. 

The integrated intensities of the 2.6 mm line of 
CO, Wen , have been derived from the composit e 
survey by iDame. Hartmann. fc Thaddeus! f|200lh . 
The data have been filtered with the moment- 
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masking technique in order to reduce the noise 
while keeping the resolution of the original data. 

Figure Q] shows the velocity-longitude profile of 
Hi emission in our region of interest (ROI). The 
preparation of maps accounting for the different 
Galactic structures present along the line of sight 
is similar to that described in detail inlAbdo et al.l 



(|2010af ) and based on a sequence of three steps: 



1. preliminary separation within Galactocen- 
tric rings; 

2. adjustment of the boundaries based on the 
velocity structures of the interstellar com- 
plexes; 

3. correction for the spill-over due to the veloc- 
ity dispersion of the broad H i lines between 
adjacent regions. 

Four regions were defined in Galactocentric dis- 
tance, namely the Local arm (Galactocentric ra- 
dius R < 10 kpc), the interarm region (R = 
10-12.5 kpc), the Perseus arm (i? = 12.5-16 kpc) 
and the region beyond the Perseus arm (which 
hosts a faint segment of the outer arm; R > 
16 kpc). The boundaries separating these regions 
under the assu mption of a flat rotation curve 
(IClemendll985h for the case of Rq = 8.5 kpc and 
6q = 220 km s _1 (where Rq and 9q are the Galac- 
tocentric radius and the orbital velocity of the lo- 
cal group of stars, respectively) are overlaid in Fig- 
ure Q] 

The preparation of the Hi and CO gas maps 
started from these preliminary velocity bound- 
aries, which were then adjusted for each line of 
sight to the closest minimum in the H i spectrum 
0. Then, the spill-over from one velocity interval 
to the next ones due to the velocity dispersion for 
the broad H i lines was corrected by fitting each H i 
spectrum with a combination of Gaussian profiles. 
We believe that this separation procedure provides 
more accurate estimates of the ISM column den- 
sities of each Galactic region than a simple slicing 
based on the rotation curve. 

In particular, effort was put into separating 
the outer arm structures from the more massive 
Perseus arm component, especially at I > 235° 



where the H i lines from the two regions merge into 
a single broad component. For directions where a 
minimum in the H i brightness temperature pro- 
file was not found near the R = 16 kpc velocity 
boundary, we integrated the profiles on both sides 
of the R = 16 kpc velocity boundary to estimate 
the Perseus and outer arm contributions. Then we 
inserted a line in the Gaussian fitting at the outer- 
arm velocity extrapolated from the I — v trend ob- 
served at I > 235° to correct for the spill-over 
due to the velocity dispersion. Given these dif- 
ficulties we expect large systematic uncertainties 
in the outer-arm 7V(H i) column densities and the 
corresponding 7-ray emissivitics will not be con- 
sidered for the scientific interpretation. We note 
that the impact on the emissivitics associated with 
the inner regions is small, < 10% as described in 

The resulting maps are shown in Figure[2]and[3] 
They exhibit a low level of spatial degeneracy, and 
thus allow us to separate the 7-radiation arising 
from the interaction with CRs in each component. 

2.1.2. Interstellar reddening 

It has been long debated whether the combi- 
nation of H 1 and CO surveys traces total column 
densities of neutral interstellar matter. By com- 
paring gas line surveys, the 7-ray o bservations by 



EGRET and dust thermal emission, iGrenier et al 



(|2005l ) reported a considerable amount of neutral 
gas at the interface between the two H 1 and CO 
emitting phases, associated with cold dust but not 
properly traced by H 1 and CO observations. Their 
finding was then confirmed by LAT data for the 



Gould Belt in the second quadrant (jAbdo et al 
2010ah . 

In order to complement the H 1 and CO 
maps, we have prepared a map derived from th e 
E(B — V) reddening map bv lSchleeel et~all (Il998h . 
The residual point sources at low latitudes were 
masked by setting to zero regions of 0.2° radius 
centered on the positions of potential IRAS point 
source H if the E(B - V) magnitude exceeded by 
> 20% that in surrounding pixels. The masked 
regions were t hen restored th rough an inpaint- 
ing technique (jElad et al.l 12005). In the course of 
the work, various source masking techniques have 



1 The minima are unlikely to be due to self absorption, 
because the velocity-distance relation is single valued in 
the outer Galaxy. 



2 http : //cdsarc . u-str asbg . f r/viz-bin/Cat?II/274 Sec 
Bcic hman et al.l lll988f) . 
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been used with negligible impact on the H i and 
CO emissivity results. 

The resulting map was fitted with a linear com- 
bination of the set of N(K i) and Wqo maps de- 
scribed in § 2.1.1. The operation was repeated for 
different choices of H i spin temperature. The fit 
was performed over the same region as for the 7- 
ray analysis, excluding a 3° x 3° region centered 
aroun d Canis Major OB 1 ( Mel'nik fe Efremovl 
199 5) where the temp erature correction applied by 
Schleeel et all (|l998l ) to construct the E(B - V) 
map from the dust thermal emission is highly un- 
certain. A preliminary fit had led to extremely 
negative residuals (< —1 mag) around I = 224°, 
b = —3°. Therefore, the residual E(B — V) map 
was calculated masking this region in the fit. We 
ar e aware that the tem perature corrections used 
by ISchlegel et al. ( 1998 ) are less reliable with de- 
creasing latitude, but the improvement we find in 
the 7-ray fit by adding the dust residual map sup- 
ports the use of their map at low latitude. 

The residual E[B — V) Tes map, after subtract- 
ing the linear combination of iV(Hi) and Wqo 
maps, is shown in Figured] (left panel). The resid- 
uals typically range from —0.5 to +0.5 magnitude. 
Large regions of positive residuals are found along 
the Galactic plane, in association with molecu- 
lar/atomic clouds. They are expected to trace 
gas not correctly traced by Hi and CO surveys. 
A remarkable region of positive residuals is de- 
tected at intermediate latitudes around I — 245°, 
b = +17°, in a region not covered by CO sur- 
veys. It corresponds to positive residuals also 
in 7-rays (§ 13. 2p and may be due to a missing, 
but possibly CO-bright molecular cloud (already 
suggested bv lDame. Hartmann. fc Thaddeusl 2001 



discussing the completeness of their survey, see 
Figure 8 of their paper). The negative residuals 
are generally small and may result from limita- 
tions in the gas column density derivation and/or 
dust spectral variations. The dust residual map 
compares well with the 7-ray residual map ob- 
tained when using only H 1 and CO to model the 
7-ray emission (Figured! right panel). The cor- 
relation between the spatial distributions of the 
dust and 7-ray residuals is statistically confirmed 
in § 13.21 Dust and 7-rays are consistent with the 
presence of missing gas in the positive residual 
clumps. The faint "glow" of negative residuals on 
both sides of the Galactic plane is driven by the 



nearby 7V(H 1) maps and it remains even when us- 
ing the smallest possible column-densities derived 
in the optically thin case. It may suggest a small 
change in average spin temperature from the mas- 
sive, compact clouds sampled in the plane to the 
more diffuse envelopes sampled off the plane, or 
it may be due to the presence of more missing 
gas in the plane than our templates can provide 
for in the fit. The dust-to-gas ratio as well as 
the 7-ray emissivity in the H 1 components would 
then be driven to higher values by the low lati- 
tude data and would slightly overpredict the data 
off the plane. 

The interpretation of the E{B — V) rcs map in 
this region of the sky is complicated by the lack 
of distance information for the dust emission. It 
is not possible to unambiguously assign the resid- 
uals to any of the regions under study. Since we 
aim at separating different regions along the lines 
of sight to investigate the CR density gradient in 
the outer Galaxy, using the H 1 and CO lines is 
essential. We have therefore used the E{B — l^) res 
map to correct for the total gas column densities. 
This approach is supported by the correlation we 
find between the dust and 7-ray data (\ 13. We 
also note that, since the dust contribution linearly 
correlated with the Hi and CO maps has been 
removed in the E(B — V) res map, this procedure 
allows us to extract the 7-ray emissivities that are 
actually correlated with the Hi and CO compo- 
nents. 

2.2. IC and Point Sources 

To model 7-ray emission not related with inter- 
stellar gas, we referred to the GALPROP code 



Strong fc Moskalenko 1998; IStrong et al 



20071 ) for 7-rays produced through IC scattering 
and to the fir st Fermi LAT cat alog (1FGL) for 
point sources (jAbdo et al.ll2010bl ). 

GALPROFH (|Strong fc Moskalenkoll99"8l : IStrong et al 



2007) is a numerical code which solves the CR 
transport equation within the Galaxy and pre- 
dicts the 7-ray emission produced via interactions 
of CRs with interstellar matter (nucleon-nucleon 
interaction and electron Bremsstrahlung) and low- 
energy photons (IC scattering). IC emission is 
calculated from the distribution of (propagated) 
electrons and the interstellar radiation fields de- 



' http://galprop.stanford.edu 
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vcloped by iPorter et al.l (|2008). Here we adopt 
the IC model map produced in the GALPROP 
run 54_77Xvarh7S in which the CR electron spec- 
trum is a djusted to agree wit h that measured by 



the LAT (jAbdo et al.H2009bh . This GALPROP 



model has been used in publications by the LAT 
collaboration such as lAbdo et al.l (|2010cJ) . 

The 1FGL Catalog is based on the first 11 
months of the science phase of the mission and 
contains 1451 sources detected at a significance 
> 4cr (the threshold is 25 in term of test statistic, 
TS0 )■ For our analysis we considered 21 point 
sources in the ROI with TS larger than 50. 

2.3. Gamma-Ray Analysis Model 

Following a well-established a pproach that 



dates back to the COS-B era (e.g., iLebrun et al 



19831 ). we modeled the 7-ray emission as a lin- 
ear combination of maps tracing the column den- 
sity of the interstellar medium. This approach 
is based on a simple, but very plausible assump- 
tion: 7-rays are generated through interactions of 
CRs and the interstellar gas, and the ISM itself is 
transparent to 7-rays. Then, assuming that CR 
densities do not significantly vary over the scale 
of the interstellar complexes under study and that 
CRs penetrate clouds uniformly to their cores we 
can model the 7-ray intensities to first order as 
a linear combination of contributions from CR 
interactions with the different gas phases in the 
various regions along each line of sight. 

We also added the IC model map by GAL- 
PROP and models for point sources taken from 
the 1FGL Catalog as described in § 2.2. To rep- 
resent the extragalactic diffuse emission and the 
residual instrumental background from misclassi- 
fied CR interactions in the LAT detector, we also 
added an isotropic component. CR interactions 
with ionized gas are not explicitly included in the 



The test statistic is defined as 



TS = 2(lnX - ln-Lo) 

where L and Lq is the maximum likelihood with and with- 
out including the source in the model, respectively. L is 
conventionally calculated as ln(L) = E;niln((?;) — T,i&i, 
where n^ and 0; are the data and the model-predicted 
counts in each pixel denoted by t he subscript i, respec- 
tively (See, e.g.. lMattox et al.|[l996h . TS is expected to be 
distributed as a \ 2 with n — no degrees of freedom if the 
numbers of free parameters in the model are respectively n 
and no (4 for sources in the 1FGL Catalog). 



model. The mass column densities of ionized gas 
are poorly known, but their contribution is gen- 
erally lower (< 10%) than that of the neutral gas 
and its scale height is much larger (~ 1 kpc com- 
pared with - 0.2 kpc; ICordes fe Lazloll2002l ). We 
therefore expect the diffuse 7-ray emission origi- 
nating from ionized gas to be largely accommo- 
dated in the fit by other components with large 
angular scales, such as the isotropic and IC ones, 
and to minimally impact the determination of the 
neutral gas emissivitics. 

Therefore, the 7-ray intensities 1^(1, b) (s _1 cm~ ; 
can be modeled as: 



sr- 1 MeV- 1 ) 



I^l,b) = 

q m>i ■ iV(Hi)(Z, b) t + £•=! q CO ,i ■ W CO {I, &)< 
+«EBv-E(B-V) res (Z ! b) + iicO,&) 



Et 



(i) 



where sum over i represents the combination of 
the Galactic regions, qm.i (s -1 sr -1 MeV -1 ) and 
<7ccm ( s_1 cm ~ 2 sr_1 MeV" 1 (K km s -1 ) -1 ) are 
the emissivitics per Hi atom and per Wco unit, 
respectively, qebv (s _1 enr 2 sr -1 MeV -1 mag -1 ) 
is the emissivity per unit of the E(B — V) rcs map 
(for which independent normalizations are allowed 
between the positive and negative residuals; see 
§ 13. 2[) . I\c and I\ S o are the IC model and isotropic 
background intensities (s -1 enr 2 sr -1 MeV -1 ), 
respectively, and PSj represents the point source 
contributions. Co mpared to the EGRET study by 
Digel et alJ (|200lh . we use two additional maps to 
better trace the ISM: the CO map in the Perseus 
arm and the E(B — V) rcs map. 

3. Data Analysis 

3.1. Observations and Data Selection 

The LAT on board the Fermi Gamma-ray 
Space Telescope, launched on 2008 June 11, is 
a pair-tracking telescope, detecting photons from 
- 20 MeV to more than 300 GcV. Details on the 
LAT instrument and prc-launch e xpectations of 



the p erformance can be found in lAtwood et al 
(12009), and the on-orbit calibration is described 



in lAbdo et all (|2009d ). 

Routine science operations with the LAT 
started on 2008 August 4. We have accumulated 
events from 2008 August 4 to 2010 February 4 to 



G 



study diffuse 7-rays in our ROL During this time 
intervai the LAT was operated in sky survey mode 
nearly all of the time, obtaining complete sky cov- 
erage every two orbits and relatively uniform expo- 
sures over time. We used the standard LAT analy- 
sis software, the Science Tools and selected events 
satisfying the standard low-backgr ound event se 



lectio n (the so-called Diffuse class; lAtwood et al 



20091 ) pi We also required the reconstructed zenith 
angles of the arrival direction of photons to be less 
than 105° and the center of the LAT field of view 
to be within 52° from the zenith, in order to re- 
duce the contamination of photons from the Earth 
limb. In addition, we excluded the period of time 
during which t he LAT detected bri ght GRBs, i.e., 
GRB080916 C (lAbdo et al.1 l2009dh . GRB090510 
( Abdo et al.1 l2009eh. GR B090902B (lAbdo et al 



2009J), and GRB090926A (|Abdo et al.ll2010eh 



3.2. Analysis Procedure 

The model described by Equation (1) was fit- 
ted to the data using the Science Tools, which 
take into account the energy-dependent instru- 
ment point spread function and effective area. We 
have analyzed the LAT data from 100 MeV to 
25.6 GeV using 13 logarithmically-spaced energy 
bands from 100 MeV to 9.05 GeV, and a single 
band above 9.05 GeV. We then have compared 
the model and data in each energy band using a 
binned maximum-likelihood method with Poisson 
statistics (in 0.25° x 0.25° bins); we thus did not 
assume an a-priori spectral shape of each model 
component except for the IC emission. For the 
other components the convolution with the instru- 
ment response functions was performed assuming 
an E~ 2 spectrum, and the integrated intensities 
were allowed to vary in each narrow energy bin. 
Changing the fixed spectral shape index over the 
range from —1.5 to —3.0 has a negligible effect 
on the obtained spectrum. In the highest energy 
band, we have set both the normalization and the 
spectral index free to accommodate the wider bin 
width. We used a post-launch response function, 
P6_V3.DIFFUSE, developed to account for the 7- 
ray detection inefficiencies due to pile-up and acci- 



Data and software are publicly avail- 
able from the Fermi Science Support Center 
(http://fermi.gsfc.nasa.gov/ssc/). For this analy- 
sis we used the P6 Diffuse selection and the Science Tools 
version v9r!6p0. 



dental coincidence in the LAT (jRando et al.ll2009l) . 
We stopped at 25.6 GeV since the photon statis- 
tics does not allow us to reliably separate different 
gas components above this energy. 

We started with point sources detected with 
high significance (TS > 100) in the 1FGL Cat- 
alog; we have 14 sources in our ROI for which 
the normalizations are set free. We also included 
8 sources lying just outside (< 5°) of the region 
boundaries, with all the spectral parameters fixed 
to those in the 1FGL Catalog. As a starting point 
we used Hi maps prepared for Ts = 125 K. We 
added model components step by step as described 
below. 

We first fitted the LAT data using Equation (1) 
without the E{B — V) rcs map and the CO map in 
the Perseus arm, and then included the CO map 
and confirmed that the fit improved significantly; 
the test statistic summed over 14 bands with sep- 
arate fits in each band (i.e., 14 more free param- 
eters) is 187.6. The 7-ray emission associated to 
the gas traced by CO in the Perseus arm is thus 
significantly detected by the LAT. 

Next we included the E(B — V) ICS map in the 
analysis. We allowed the independent normaliza- 
tions between the positive part and the negative 
part of the E{B — V) ICS map, and found that the 
normalizations differ with each other. We thus will 
use the independent normalizations hereafter. We 
chose this model to better represent the LAT data 
and constrain the CR distributions, and leave a 
detailed discussion about the use of dust as ISM 
tracer to a dedicated paper. The improvement 
of the fit is very significant: TS = 1119.6 for 28 
more free parameters. The correlation between 
the E[B — V) rcs map and the 7-ray residual map 
obtained by the fit without the E{B — V^res map, 
shown in Figure SI further supports the use of 
E(B — V) Tes map in our analysis. 

We also tried a fit without the IC component 
to assess the systematics. The effects on the de- 
rived emissivities are typically 2-3% and ~ 5% 
for qui and qco, respectively. They are much 
smaller than the statistical errors and systematic 
uncertainties (see below), although the inclusion 
of the IC map improves the fit to the LAT data. 
Therefore the uncertainties on the IC model have 
no significant impact on our analysis due to its 
rather flat distribution across the region of interest 
while the gas in the ISM is highly structured. On 
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the other hand, lowering the threshold for point 
sources down to TS = 50 yields an about twice 
smaller emissivity for the Wco map in the Perseus 
arm. The emissivities of other components are un- 
changed within the statistical errors. This is plau- 
sibly due to the very clumpy distribution of the 
clouds in the Perseus arm as seen by a terrestrial 
observer, see Figure which makes it difficult to 
separate from that of some discrete sources. We 
thus use Equation (f ) with point sources detected 
at TS > 50 in the 1FGL Catalog H as our baseline 
model, but wc do not consider the highly uncer- 
tain CO emissivities in the Perseus arm for the 
discussion. 

We summarize the results in Tables [T] and [5] for 
Tg = f 25 K and 250 K, respectively, and the num- 
ber of counts in each energy bin in Table [3J The 
differential emissivities are multiplied by E 2 where 
E is the center of each energy bin in logarithmic 
scale. They are given for each model component. 
We note that our isotropic term (ii so ) includes the 
contribution of the instrumental background and 
might partially account also for ionized gas (see 
§ 12 . 3[) . thus it is significantly larger tha n the extra 



galactic diffuse emission reported by lAbdo et al 
(<20f0dl) . 



To illustrate the fit quality, we give the data 
and model count maps and the residual map in 
Figure [5] (for Tg = 125 K), in which residu- 
als (data minus model) are expressed in approx- 
imate standard deviation units (square root of 
model-predicted counts). Although some struc- 
tures (clustering of positive or negative residuals) 
are observed, the map shows no excesses below 
—4 a and above 6 a. Over 99% of the pixels are 
within ±3 a. We thus conclude that our model 
reasonably reproduces the data. 

Figure [6] presents the fitted spectra for each 
component. The emission from the Hi gas domi- 
nates the 7-ray flux. Although the emission from 
the gas in the CO-bright phase and that traced by 
E(B — V) rcs is fainter than the IC and isotropic 
components, their characteristic spatial structures 
(see Figure [2] and [3]) allow their spectra to be reli- 
ably constrained. 

To examine the effect of the optical depth cor- 



rection applied to derive the Hi maps, as antici- 
pated above we tried several choices of a uniform 
Tg. We stress that the true Ts is likely to vary 
within clouds, but we stick to this simple ap- 
proximation exploring the following values: 100 K 
(which is a reasonable lower limit in the uniform 
approximation^, 250 K and 400 K (which are the 
two values indicated by absorption measur ements 
in the outer Galaxy by IDickev et al. l l2009lM and 
the optically-thin approximation (which yields 
the lower limit allowed on the atomic column 
densities). The results on the maximum log- 
likelihood values are summarized in Table |4] to- 
gether with the integrated H i emissivities ob- 
tained above 100 MeV in each region. The evolu- 
tion of ln(L) with Tg is plotted in Figure [7J The 
Hi emissivity varies by +15/ — 10% for the Lo- 
cal arm, +10/ — 0% for the interarm region, and 
+50/ — 25% for the Perseus arm with respect to 
the Tg = 125 K case. We observe an increase of 
ln(L) with increasing spin temperature. Consid- 
ering the fact that Tg ~ 250 K is a typical value 
in the second quadrant of t he outer Galaxy ac - 
cording to a recent study bv IDickev et al. ( 20091 ) 
and because ln(L) saturates at Tg > 250 K, we 
regard 250 K as a plausible estimate of the aver- 
ag e Tg in our ROI- U nfortunately, the estimates 
bv IDickev etldl (|2009h have a rather large uncer- 
tainty (about ±50 K) in each Galactocentric ra- 
dius bin, and they do not cover the region in the 
th ird quadrant we ar e investigating (sec Figure 5 
of IDickev et~aTll2009h . In the following sections, 
we will concentrate on Tg = 125 K for compar- 
ison with previous 7-ray measurements and on 
Tg = 250 K which agrees well with H 1 absorption 
and the LAT data. 

4. Results and Discussions 

4.1. Emissivity Spectra of Atomic Gas 

In Figures [8] and [9] (left panels) , wc report the 
emissivity spectra found per H atom in the Local 
arm, interarm, Perseus arm and outermost regions 
for Tg = 125 and 250 K, respectively. For compar- 
ison with the local interstellar spectr um (LIS) we 
also plot the model spectrum used in lAbdo et al.l 



6 Spectral parameters of point sources of TS=50-100 are 
fixed to those given in the 1 FGL Catalog in the highest 
energy bin. (9.05-25.6 GeV) 



7 A truncation at 95 K was applied for channels where the 
brightness temperature was larger. 

8 Note that, however, the data used bv IDickev et al. 
do not cover the third Galactic quadrant 
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(|2009al ) which agrees well with LAT data in a mid- 
to-high-latitude region (22° < |6| < 60°) of the 
third quadrant (assuming Ts = 125 K). We see 
that the spectral shape of the Local-arm emissiv- 
ity agrees well with the model for the LIS and does 
not depend on the choice of spin temperature. The 
integral emissivity of the Local arm is 10 % lower 
than that reported bv lAbdo et all (|2009al ) for the 
same spin temperature. This difference is not sig- 
nificant given the uncertainties in the kinematic 
separation of the gas components. The present re- 
sult is also consist ent with the measur ement in the 
second quadrant ( Abdo et al.l l2010ah . Together 
they show that the CR density along the Local 
arm is rather uniform within 1 kpc around the 
Sun, both in the second and third quadrants. 

The comparison of the data with the model 
emissivity expected for the Local arm region based 
on locally-measured CRs (Figures E] and indi- 
cates a better fit for higher Ts; Ts = 125 K gives 
emissivities 15-20% lower than the model, whereas 
Ts = 250 K shows better agreement by about 10%. 
Although the theoretical emissivity has uncertain- 
ties due t o imperfect knowl edge of the CR spec- 
trum fsee lAbdo et al.ll2009al ). the fact that a high 
Ts value yields a better match both to the local 
absolute emissivity and to the spatial distribution 
of the diffuse emission (Figure [7|) leads to larger 
Ts than a value conventionally used in 7-ray as- 
trophysics (125 K). This is in accord with inde- 
pendent estimates of Ts as discussed in § 13.21 

We also observe that the emissivity spectra 
do not vary significantly with Galactoccntric dis- 
tances in the outer Galaxy. To examine the spec- 
tral shape more quantitatively, we present the 
emissivity ratios of the interarm and Perseus re- 
gions relative to the Local arm in the right panels 
of Figures [8] and [H The spectral shape in the in- 
terarm region is found to be consistent with that 
in the Local arm; a fit to the data for Ts = 125 K 
with a constant ratio gives \ 2 = 7.3 for 13 degrees 
of freedom. Although the fit is not fully accept- 
able for the Perseus arm (x 2 = 24.3), the large x 2 
is driven solely by the last bin. We note a possible 
interplay between the Perseus-arm and the adja- 
cent outer-arm emissivities in the highest energy 
bins (see left panels of Figure [5] and O . It can be 
due to a small but non-negligible spatial difference 
between the modeled templates and data and/or 
to the presence of unresolved point sources (gener- 



ally harder than diffuse emission) . Photon fluctua- 
tions from the structured gas components can also 
lead the fit to a slightly different solution in the 
spatial separation of the components. One would 
expect these possible systematic uncertainties to 
become important at high energy given the limited 
counts in the overall map. It is difficult to quan- 
titatively test these effects without knowledge of 
the true model distributions, but we can note that 
the small deviations seen at 400-560 MeV and 1.6- 
2.2 GeV from a constant ratio are not confirmed 
by the general trend of the other points. They in- 
dicate that there are systematic uncertainties not 
fully accounted for by the statistical errors in the 
fit. We thus do not claim nor deny the spectral 
softening of the Perseus arm at high energy. A 
test using Ts = 250 K for the iV(Hi) maps gives 
the same conclusion on the spectral shape. We 
thus conclude that the spectral shapes are consis- 
tent with the LIS in the 0.1 — 6 GeV energy band, 
independent of the assumed Tg. Considering that 
these 7-rays trace CR nuclei of energies from a few 
GeV to about 100 GeV (see, e.g., Figure[n]of[Morj| 
19971 ). LAT data indicate that the energy distri- 
bution of the main component of Galactic CRs 
does not vary significantly in the o uter Galaxy 
i n the third quadrant. We note that lAbdo et al 



(|2010al ) reported a possible spectral hardening in 
the Perseus arm in the second quadrant. This 
might be due to the presence of the very active 
star-forming region of NGC 7538 and of CRs hav- 
ing not diffused far from their sources, or to con- 
tam ination by hard un resolved point sources. In 
fact I Abdo et al. ( 2010a ) did not rule out the possi- 
bility that their result is due to systematic effects. 

4.2. Calibration of Molecular Masses 

High-energy 7-rays are a powerful probe to de- 
termine the CO-to-H2 calibration ratio, Xco, if 
the CR flux is comparable in the different gas 
phases inside a cloud. Since the 7-ray emission 
from the molecular gas is primarily due to CR in- 
teractions with H2, and since the molecular bind- 
ing energy is negligible in processes leading to 7- 
ray production, the emissivity per H2 molecule is 
expected to be twice the emissivity per H 1 atom. 
Then, under the hypothesis that the same CR 
flux penetrates the Hi- and CO-bright phases of 
an interstellar complex, we can calculate Xqo as 
qco = 2Xco • 9hi- 
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We show gco as a function of qm for the Lo- 
cal arm and the interarm region in Figure 1101 
We do not consider the correlation in the Perseus 
arm, because qco from this region is affected by 
large systematic uncertainties (see § 3.2). Since 
the emissivity associated with the CO-bright gas 
is not well determined in the lowest energy range 
(see Table Q] and [2} because of the poor angular 
resolution of the LAT, and the fit at very high en- 
ergy is affected by larger uncertainties (§ 4.1), we 
have plotted only data in the 200 McV-9.05 GeV 
range. The linear correlation supports the as- 
sumption that Galactic CRs penetrate molecular 
clouds uniformly to their cores. It also suggests 
that contamination from point sources and CR 
spectral variations within the clouds are small. 

We have derived the maximum-likelihood esti- 
mates of the slope and intercept of the linear re- 
lation between qco and qm taking into account 
that qco an d qm are both measured (not true) val- 
ues with known uncertainties. The resulting inter- 
cepts are compatible with zero. The Xco values 
we have obtained for T s = 250 K are (2.08±0.11) x 
10 20 cm -2 (K km s" 1 )" 1 for the Local arm (R < 
10 kpc) and (1.93± 0.16) x 10 20 cm- 2 (K km/s)" 1 
for the interarm region (R = 10-12.5 kpc). De- 
creasing the spin temperature to 125 K does not 
affect the Xco derivation in the well resolved, not 
too massive, clouds of the Local arm where we find 
X CO = (2.03 ± 0.11) x 10 20 cm" 2 (K km s" 1 )" 1 . 
On the other hand, the separation in the 7-ray 
fit between the dense H 1 peaks and clumpy CO 
cores becomes more difficult for more distant, 
less resolved clouds where Hi and CO tend to 
peak in the same directions. A change in the 
largest iV(Hi) column-densities from the optical 
depth correction can impact the Xco determina- 
tion in two ways: first by changing the qui emis- 
sivity, second by modifying the N(H 1) contrast 
within the cloud, hence the Hi and CO separa- 
tion. The global impact is mild since we find 
Xco = (1-56 ± 0.17) x 10 20 cm- 2 (K km s" 1 )" 1 
in the interarm region for Tg = 125 K. 

Abdo et al. ( 2010al ) reported comparable val- 
ues of Xco in the second quadrant for Tg = 125 K: 
they obtained (1.59±0.17) x 10 20 cuT 2 (K km s" 1 )" 
and (1.9 ± 0.2) x 10 20 cur 2 (K km/s)" 1 for the 
Local arm (R < 10 kpc) and the Perseus arm. 
Given the systematic uncertainty in Xco ; roughly 
of the order of 0.3 x 10 20 cur 2 (K km s -1 ) -1 , 



due to H 1 optical depth correction, the results of 
both studies point to a rather uniform ratio over 
several kpc in the outer Galaxy. Yet, these values 
are twice larger than found in the very nearby 
Gould-Belt clouds of Cassiopeia and Cepheus, 
Xco = (0.87 ± 0.05) x 10 20 cm- 2 (K km/s)- 1 , 
Tg = 125 K. However, we confirm that the in- 
crease in Xco beyond the solar circle is signifi- 
ca ntly lower than the tr end adopted in the model 
of IStrong et all (|2004bl ). What fraction of the 
Gould-Belt to Local-arm difference in the average 
Xco can be attributed to a difference in the spa- 
tial sampling (resolution) of the clouds remains to 
be investigated. 

Nearly the same region has been analyzed 
by iDigel et al.l (|200ll ) using EGRET data. The 
main difference from their analysis is our im- 
proved scheme for the kinematical separation 
of the ISM components along the lines of sight 
and the inclusion of the reddening residual map. 
The Xco value measured in the Local arm by 
EGRET, (1.64 ± 0.31) x 10 20 cm- 2 (K km s" 1 )" 1 
(for Tg = 125 K), is statistically compatible with 
ours. The fact that we excluded the region of 
Monoccros R2 from the analysis can also explain 
in part this difference. 

Because of the pile-up of different clouds along 
a line of sight, the derivation of individual cloud 
masses is beyond the scope of this study. Let us 
just note that Maddalena's cloud, with its very 
low rate of star formation, seems to share a quite 
conventional Xco factor. Further investigation, 
including higher resolution 7-ray maps when more 
high-energy LAT data become available, is re- 
quired to fully understand the mass distribution 
in the clouds. 

4.3. The Gradient of CR Densities beyond 
the Solar Circle 

In Figure [TT] (left panel) we show the emissivity 
gradient found beyond the Solar circle for different 
spin temperatures. Here we do not include the re- 
sults for the optically-thin approximation which is 
equivalent to an infinitely high Tg and gives similar 
emissivities to Tg = 400 K. The typical statistical 
errors associated with these measurements are il- 
lustrated in the right panel for the Tg = 125 K 
case. In the right panel, a shaded area shows 
the characteristic systematic error due to the LAT 
event selection efficiency, evaluated to be ~ 10% 
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in the energy range under study. 

In order to evaluate the impact of the deli- 
cate separation of the gas in the outermost re- 
gion, we have compared two extreme cases. The 
first one adopts the kinematic R — 16 kpc bound- 
ary and applies no correction for velocity disper- 
sion, the second assigns all the outer-arm gas to 
the Perseus arm. The emissivity in the Perseus 
arm differs by about 5% from the original one, 
and those in the Local arm and interarm regions 
hardly change. Therefore these effects are signif- 
icantly smaller than the uncertainties due to the 
optical depth correction of the H i data. We also 
note that the main effect of the LAT selection ef- 
ficiency uncertainty is to rigidly shift the profile 
without any significant impact on the gradient. 

We thus conclude that the most important 
source of uncertainty in the CR density gradient 
derivation is currently that in the N(H i) determi- 
nation. This is mainly because the optical depth 
correction is larger for dense H i clouds in the Lo- 
cal and Perseus arms than for diffuse clouds in 
the interarm region. The loss in contrast between 
the dense (low-latitude) and more diffuse (mid- 
latitude) H i structures resulting from an increase 
in spin temperature affects the fit, particularly in 
the Perseus component which is more narrowly 
concentrated near the plane. When probing the 
CR densities as the "ratio" between the observed 
numbers of 7-rays to H atoms, at the precision 
provided by the LAT the uncertainties in the ISM 
densities arc dominant. 

4-3.1. Comparison with EGRET and the arm /interarm 
contrast 

An in teresting finding of the former EGRET 
analysis (jDigel et al 1 l200lh was an enhancement of 
the 7-ray emissivity in the Perseus arm compared 
with the interarm region. This possibility is rel- 
evant for models of diffuse 7-ray emissions based 
on the assump tion that CR and IS M densities are 
coupled (e.g.. Irlunter et al.lll997l and references 
therein). 

The Local arm emissivity obtained by the 
EGRET study for T s = 125 K is (1.81 ± 
0.17) x 10 -26 photons s _1 sr _1 H-atom -1 , which 
is ~ 25% larger than our LAT result. However, 
the two studies are based on different H 1 surveys 
which yield different total N(H 1) column densities 



integrated along the lines of sight. The column- 
density ratios between the surveys varies from 0.6 
to 1.0 within the ROI, with an average value of 
~ 0.8. The difference is likely due to the improved 
correction for stray-r adiation in the more recent 
survey, as discussed in iKalberla et all ([20051) . The 
EGRET Local-arm emissivity scaled by 0.8 is in 
good agreement with our result for the same spin 
temperature. If we do not include the E(B — V) leB 
map in the fitting, we obtain an emissivity of 
(1.68 ±0.05) x 10" 26 photons s" 1 sr" 1 H-atom" 1 
which is still consistent with the down-scaled 
EGRET result within -15%. We thus conclude 
that our result is consistent with the previous 
study but is more reliable because of higher 7-ray 
statistics, finer resolution and an improved H 1 gas 
survey. 

We can therefore compare the present emis- 
sivity gradient (for consistency in the case Tg = 
125 K) with that reported by the EGRET study, 
as summarized in Figure [TTJ (right panel in which 
the EGRET results multiplied by 0.8 are also 
shown) . Although we observe good agreement be- 
tween the two studies in the Local and the Perseus 
arms, this is not true for the arm/inter arm con- 
trast. The difference could be due to the simple 
partitioning in cloud velocity used for the EGRET 
study. The H 1 mass obtained for clouds in the 
interarm region with the simple partitioning is 
20% - 40% larger (for T s = 125 K) than with 
our separation scheme, exaggerating the amount 
of gas in the interarm region, and thus lowering 
the emissivity by the same amount. Our emis- 
sivity profile is thus consistent with the previous 
study, but with improved precision (smaller sta- 
tistical errors) and accuracy (more reliable region 
separation method and better estimation of the 
point source contributions). We thus do not con- 
firm a marked drop in the interarm region. 

Low spin temperatures yield a smooth de- 
cline in H 1 emissivity to R ~ 16 kpc in the 
outer Galaxy, without showing a significant cou- 
pling with ISM column densities. The Perseus- 
to-interarm contrast is at most of the order of 
15%-20% for high spin temperatures as shown 
in the left panel of Figure [TT] These profiles are 
similar at all energies, in particular at high en- 
ergies where the component separation is more 
reliable thanks to the better angular resolution. 
The surface density of H 1 in the Perseus arm is 
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on average 30%-40% higher than in the interarm 
region. Therefore, even if we adopt Ts = 400 K 
which gives the largest arm-inter arm contrast, the 
coupling scale (or the co upling length) betw een 
the CRs and matter (e.g., iHunter et al.lll9971 ) re- 
quired to agree with the LAT data would be larger 
than those usually a ssumed for this ty pe of model 
(~ 2 kpc. See e.g.. iDigel et~aD l200ll Figure 7). 
Whether the true emissivity profile exhibits a 
small contrast between the arms or smoothly de- 
clines with distance is beyond our measurement 
capability without further constraints on the H 1 
column-density derivation. New H 1 absorption 
measurements will allow us to investigate this is- 
sue with better accuracy. 

4-3.2. Comparison with a Propagation Model: 
the CR gradient problem 

To compare wit h the second quadrant results 
(jAbdo et al.ll2010al ). we have integrated the emis- 
sivities above 200 MeV for T s = 125 K. We find 
values of (0.817±0.016) x 10" 26 photons s" 1 sr" 1 H-a 
(0.705±0.018)xl0" 26 photons s" 1 sr" 1 H-atom" 1 
and (0.643±0.022)xl0~ 26 photons s" 1 sr" 1 H-atom" 
for the Local-arm, interarm, and Perseus-arm re- 
gions, respectively. The nearer value is about 20% 
lower than in the 2nd quadrant (which, however, 
samples very nearby clouds in the Gould Belt) and 
the outer ones compare very well with the 2nd 
quadrant measurement over the same Galactocen- 
tric distance range. Despite the uncertainties due 
to the optical-depth correction (that might have a 
different impact in the two quadrants), both LAT 
studies consistently point to a slowly decreasing 
emissivity profile beyond R = 10 kpc. 

Let us consider the predictions by a CR prop- 
agation model to see the impact of such a flat 
profile on the CR source distribution and prop- 
agation parameters. We adopted a GALPROP 
model, starting from the configuration used for 
the run 54_77Xvarh7S which we used to predict 
the IC contribution. The CR source distribution 
in this model is, 

where Rq — 8.5 kpc is the distance of the Sun 
to the Galactic center. As shown in Figure [T2| 
this function is intermediate between the distribu- 
tion of supernova remnants (SNRs) obtained from 



the Yi-D relation (jCase fc Bhattacharva [1998) 
and that traced by pulsars (jLorimerl 120041 ) . The 
boundaries of the propagation region are cho- 
sen to be i?h = 30 kpc (maximum Galactocen- 
tric radius) and Zh = 4 kpc (maximum height 
from the Galactic plane), beyond which free es- 
cape is assumed. The spatial diffusion coeffi- 
cient is assumed to be uniform across the Galaxy 
and is taken as D xx = (3Do(p/4: GV) 5 , where 
j3 = v/c is the velocity of the particle relative 
to the speed of light and p is the rigidity. We 
adopted D = 5.8 x 10 28 cm 2 s" 1 and 5 = 0.33 
(Kolmogorov spectrum). Rcacceleration due to 
the interstellar magnctohydrodynamic turbulence, 
which is thought to reproduce the observed B/C 
ratio at low energy, assumes an Alfven velocity 
V\ — 30 km s _1 . The CR source distribution and 
propagation model paramete rs have been used of - 
ten in the literature (see e.g., Strong et al. 2004al ). 
We note that the same CR source distribution and 
similar propagation para meters are adopted in the 
GALPROP run used bv lAbdo et ail (|2010af) . 

The left panel of Figure [13] compares the cal- 
culated profile (solid line) with LAT constraints 
(bow-tie plot bracketing the profiles obtained for 
different Ts ; see the left panel of Figure HT|) . The 
model is normalized to the LAT measurement in 
the innermost region. Despite the large uncertain- 
ties, LAT data lead to a significantly flatter pro- 
file than predicted by our model; the LAT results 
indicate to a factor of 2 larger emissivity (CR en- 
ergy density) in the Perseus arm even if we assume 
T s = 100 K. The higher T s makes the discrep- 
ancy larger, hence the conclusion is robust. Not 
using the E(B — V) res map in the analysis does 
not change the conclusion, since the emissivitics 
in the interarm region and the Perseus arm are 
almost unaffected by its presence. 

The discrepancy between the 7-ray emissivity 
gradient in the Galaxy and the distribution of pu- 
tative CR sources has been known as t he "gradi- 
ent p roblem" since the COS-B era (e.g.. iBloemen 
1989). It has led to a number of possible interpre- 



tations, including, for the specific case of the outer 
Galaxy, the possibility of a very steep gradient i n 
Xco beyond the solar circle (jStrong et al.ll2004bl ) . 
The cmissivities in the outer Galaxy were more dif- 
ficult to determine in the COS-B/EGRET era due 
to lower statistics and higher backgrounds. Now, 
thanks to the high quality of the LAT data and the 
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improved component separation technique applied 
to gas line data, we measure a flat H i emissivity 
gradient in the outer Galaxy together with a flat 
evolution of Xco over several kpc, so the gradient 
problem requires another explanation. 

The most straightforward possi bility is a larger 
halo s iz e (zh), as discus sed b y. e.g.JStecker fc Jones 
19771) . iBloemen] (|l989h and lStrong fc MoskalenkoT 
19981 ). We therefore tried several choices of 
Zh and Dq as summarized in the dotted lines 
in the same figure. The values of Do are cho- 
sen to reasonably reproduce the LIS of protons 
and electrons, B/C ratio and 10 Be/ 9 Be ratio at 
th e solar system, and are simi lar to those given 



The exploration could be extended to other re- 
gions of the parameter sp ace or to a non-un iform 



Strong fc M oskalenko (|l998l) . All models are 



normalized to the LAT data in the Local arm. 
Models with Zh — 4 kpc or smaller arc found to 
give too steep emissivity gradients. A CR source 
distribution as in equation 2 with a very large 
halo (zh > 10 kpc) provides a gradient compati- 
ble with the 7-ray data, if we fully take into ac- 
count the systematic uncertainties. We note that 
Zh = 10 kpc is sti ll compatible with 10 Be/ 9 B e 



measurements (e.g.. IStrong fc Moskalenkdll998f ) 



Considering the large statistical and systematic 
uncertainties in the SNR distribution, a flatter 
CR source distribution in the outer Galaxy also 
could be possible. We thus tried a modified CR 
source distribution, in which the distribution is 
the same as equation 2 below i?bk and constant 
beyond it (see a thin solid line of Figure rjS] as 
an example). Figure [T3l right shows the models 
with several choices of i?bk for Zh = 4 kpc and 
D = 5.8 x 10 28 cm 2 s _1 . We obtained a reason- 
able fit to the data using a flat CR source distri- 
bution beyond R = 10 kpc. Such a constant CR 
source density in the outer Galaxy is in contrast 
not only with the (highly uncertain) distribution 
of SNRs, but also with other tracers of massive 
star formation and SNRs, like, 1) CO lines which 
trace the interstellar ph ase where massive stars 
form (e.g., Ferriere 200lh . 2) OB star counts (e.g., 



Bronfman et al. 2000l ), and 3) the 26 Al line which 



is related to the injection of stel lar nucleosynthesi s 



products in the ISM by SNRs (jDiehl et al.l l2006) 



However, a very large halo size and/or a flat CR 
source distribution just beyond the solar circle 
seem to be favored by the LAT data. 

The above discussion depends on the propaga- 
tion parameters and the solution is not unique. 



Evoli et alJboOSl ). but 



diffusion coefficient (e.g. 
examining propagation models in detail is beyond 
the scope of our study. Our bottom line is that 
the analysis of LA T data presented here and by 
Abdo et al. ( 2010a[ ) consistently show that the CR 
density gradient in the outer Galaxy is flatter than 
expectations by commonly-used propagation mod- 
els. In the future, the extension to the inner part 
and the accurate determination of the gradient 
over the whole Galaxy will be key to constrain- 
ing the CR origin and transport. 

We also note that a spin temperature T$ > 
250 K, which is fav ored by recent stud ies in the 
outer Galaxy (e.g., iDickev et al. 120091) . gives a 



small arm/interarm contrast at the 10-20% level 
that is not fully compatible with the propagation 
models (including the one we adopted here) which 
predict a monotonic CR gradient. 

Even though the present analysis includes 
a dust template to account for the abundant 
missing gas present locally at the interface be- 
tween the Hi and CO-bright phases, an alter- 
native way to reconcile the flat emissivity pro- 
file and a marked decline in CR density in the 
outer Galaxy is to invoke an increase in miss- 
ing gas mass with Galactocentric distance in 
the low metalli city environments of the outer 
Galaxy (see e.g. Papadopoulos. Thi fc Viti 2002t 
Wolfire et alj 120101 ) beyond the local correction 
applied here. We note that the large masses 
of dark gas in the outer Galaxy suggested by 
Papadopoulos. Thi fc Vitil (|2002l ) (outweighing 



that of Hi by a factor of 5-15) might explain 
our results, whereas the remarkably constant dark 
gas fraction of 30% wi th mild dependence on 
metallicity suggested by IWolfire et al. ( 2010l ) is 
not sufficient to explain the large H i emissivities 
measured by the LAT beyond the solar circle. 

5. Summary and Conclusion 

We have studied the diffuse 7-ray emission in 
the third Galactic quadrant using the first 18 
months of Fermi LAT science data. Thanks to 
the excellent performance of the LAT, we have 
obtained high-quality emissivity spectra of the 
atomic and molecular gas (traced by Wco) in the 
100 MeV - 25.6 GeV energy range. 
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At the level of accuracy allowed by the LAT, 
the study of CR densities from 7-ray observations 
is now mostly limited by the understanding of the 
ISM mass tracers, notably by the uncertainties in 
the derivation of atomic gas column densities from 
H 1 surveys and by the distribution of gas not ac- 
counted for by radio and microwave line surveys. 
In spite of those uncertainties, robust conclusions 
can be drawn concerning the ISM and CRs. 

The molecular mass calibration ratio of the Lo- 
cal arm is found to be ~ 2xl0 20 cm -2 (K km s _1 )~ 
significantly larger than that for the very local 
Gould-Belt clou ds in the second Ga lactic quad- 
rant reported by lAbdo et all (j2010al ). No signifi- 
cant differences of the ratio are found between the 
Local arm and the interarm region. 

No significant variations in the CR spectra are 
found across the outer Galaxy in the region stud- 
ied, and no large contrast in emissivity is seen in 
the interarm region between the Local and Perseus 
arms (a contrast < 10-20% is allowed by data). 
The measured gradient is much flatter than predic- 
tions by a widely-used propagation model assum- 
ing that the CR source distribution largely peaks 
in the inner Galaxy. A larger halo size and/or 
a flatter CR source distribution beyond the so- 
lar circle than those usually assumed are required 
to reproduce the LAT data, while other scenar- 
ios such as a non-uniform diffusion coefficient or 
vast amounts of missing gas in the outer Galaxy 
are also possible. Reliable determinations of the 
amount of atomic hydrogen in the plane are key 
to better constraining the property of CRs in our 
Galaxy. 
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Fig. 1. — Longitude- velocity diagram of the av- 
erage intensity of the 21 cm line (in unit of K) for 
— 15° < b < 20°. Preliminary boundaries between 
the four Galactocentric annuli are also presented 
(See § 2.1.1 for details). The lowest contour cor- 
responds to 2 K and the contour interval is 3 K. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 
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Fig. 2. — Maps of 7V(Hi) (in unit of 10 20 atoms cm~ 2 ) for the Local arm (top left), interarm (top right), 
Perseus arm (bottom left), and outer arm (bottom right) regions, obtained for a spin temperature Ts = 125 K. 
The outlined area in the bottom right corner is not used in the analysis (see § 1). The maps have been 
smoothed with a Gaussian with a = 1° for display. 
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Fig. 3. — Maps of Wco (in unit of K km s _1 ) for the Local arm (top left), interarm (top right), and 
the Perseus arm (bottom left) regions. The small box in the bottom right corner indicates the area not 
considered in the analysis. The maps have been smoothed with a Gaussian of a = 0.25° for display. 
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Fig. 4. — (left) Residual E(B — V) map in unit of magnitudes, obtained by subtracting the parts linearly 
correlated with the combination of iV(Hi) and Wco maps. The small box in the bottom right corner shows 
the area not considered in the analysis. The map has been smoothed with a Gaussian of a = 0.25° for 
display, (right) 7-ray residual (data minus model) map obtained by the fit without the E(B — V) rcs map 
(only H 1 and CO maps) in unit of standard deviations (square root of model-predicted counts, saturated 
between — 3<7 and +3ct). The map has been smoothed with a Gaussian of a = 0.5°. 
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Fig. 5. — Data count map (top left), model count map (top right) and the residual (data minus model) map 
in units of standard deviations (bottom left, saturated between — 3er and +3cr) above 100 MeV obtained by 
our analysis. Point sources with TS > 50 in the 1FGL included in the fit are shown by crosses in the model 
map. Data/model count maps are in 0.25° x 0.25° pixels, and the residual map has been smoothed with a 
Gaussian of a — 0.5°. 
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Fig. 6. — 7-ray spectra over the region of interest obtained from data and from the fitted model (for each 
gas phase, IC and isotropic components, and for point sources). 
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Fig. 7. — Variations of the log-likelihood value for several choices of Ts (the scale has a fixed offset). 
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to that of the Local arm. In both panels, Ts = 125 K is assumed and spectra in the interarm region and the 
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Fig. 9.— The same as Figure U but for T s = 250 K. 
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Fig. 10. — Correlation between the Hi and CO emissivities in the 200 MeV-9.05 GeV energy range for the 
Local arm and the interarm regions. The cases of Ts = 125 K and 250 K are shown in the left panel and the 
right panel, respectively. Dotted lines show the best linear fits. Each data point corresponds to an energy 
bin used in the 7-ray analysis. (See Table Q] and [2]) 
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Fig. 11. — (left) Emissivity gradient for several choices of Ts. (right) Emissivity gradient obtained by the 
LAT compared with the EGRET results under the assumption of Ts = 125 K. The shaded area indicates the 
systematic uncertainty in the LAT selection efficiency of ~ 10%. The EGRET points have been downscaled 
by 20% to account for the change in H 1 survey data between the two studies (see § 14.3. ip . 
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Fig. 12. — The CR source distribution adopted in our ba seline GALPROP model (s olid line), compared 
with the SNR distributi on obtained by the £ — D relation (jCase fc Bhattacharvalll998l ) and that traced by 
the pulsar distribution (jLorimerl 120041 ) shown by dotted lines. The thin solid line represents an example of 
the modified distributions introduced to reproduce the emissivity gradient by the LAT. 




Fig. 13. — Comparison of the emissivity gradient obtained by the LAT and model expectations using 
GALPROP. The left panel shows models with different halo sizes and diffusion lengths: (zh,-Do) = (1 kpc, 
1.7 x 10 28 cm 2 s" 1 ), (2 kpc, 3.2 x 10 28 cm 2 s" 1 ), (4 kpc, 5.8 x 10 28 cm 2 s" 1 ), (10 kpc, 12 x 10 28 cm 2 s" 1 ) 
and (20 kpc, 18 x 10 28 cm 2 s _1 ). The solid line is for Zh = 4 kpc. The right panel shows different choices 
of the break distance beyond which a flat CR source distribution is assumed: i?bk = 10-14 kpc with 1 kpc 
steps. 
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Table 1 

a summary of fit parameters with 1 sigma statistical errors, under the assumption of 

T s = 125 K. 
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1 


13 


± 





16 





76 


± 





16 


1 


28 


± 0.06 


1.60- 


-2.26 


1 


60 


± 





08 


1 


24 


± 





08 


1 


07 


± 





Oil 


1 


9 


= 





5 


6 





± 


1 





4 


1 


± 


1 


3 


7 ± 4 





07 


± 





16 





53 


± 





16 





91 


± 0.07 


2.26- 


-3.20 


1 


26 


± 





08 


1 


02 


± 





08 





9 1 


± 





Oil 


1 


3 


± 





5 


3 


2 


± 


1 





3 


1 


± 


1 


3 


3 ± 4 





82 


± 





16 


0. 


58 


± 





16 





85 


± 0.07 


3.20 


-4.53 





80 


± 





00 





76 


± 





08 





75 


± 





09 





1 


= 







5 


7 


± 


1 





3 


5 


= 


1 


3 


11 ± 4 





7 1 


± 





10 





42 


± 





16 





93 


± 0.08 


4.53 


-6.40 





59 


± 





09 





11 


± 





08 





57 


± 





Oil 





6 


± 





5 


2 


9 


± 


1 





3 


6 


± 


1 


3 


± 





7 1 


± 





16 





02 


± 





17 





86 


± 0.08 


6.40- 


-9.05 





51 


± 





09 





31 


± 





08 





28 


± 





Oil 





9 


± 





5 


2 


2 


± 


1 





1 





= 


1 


1 


2 ± 3 





52 


± 





15 





07 


± 





16 





62 


± 0.08 


9.05- 


-25.6 





34 


± 





05 





22 


± 





05 





Oil 


± 





05 





7 


= 





3 


2 


2 


± 





6 


1 


1 


= 





.7 


3 ± 2 





07 


± 





06 





00 


± 


0. 


08 





• If) 


± 0.05 



a Units: g HIi (10~~ 24 McV^ s ~ x si" x McV" 1 ), 

1—2 ]\^,-,-\7-2 = — 1 „ m — 2 „ — 1 TV/T,-,-^— 1 man- - I'i 2 . „. fill 3 



q COi i(10~ MeV z s _ 
- cm -2 sr -1 MeV -1 ) 



q EBV (10 _z MgV z s" 1 cm~ z sr" 1 McV -1 mag' 1 ), E ■ q igo (10 -,:S MeV z s~ 

The subscripts refer to four regions defined to perform the analysis; 1) Local arm, 2) interarm region, 3) Pers 
Some parameters are not well determined and their best-fit value is consistent with 0. We present them for cc 



i 4) beyond the Perseus ; 



Table 2 

a summary of fit parameters with 1 sigma statistical errors, under the assumption of 

T s = 250 K. 



Energy 
GeV 


E 


2 


9HI 


,1 


E 


2 


'JHI 


.2 


E 


2 


'/HI 


,3 


E 2 




9HI,4 


E 2 


■ 9CO,l 


E 2 


■ 3CO,2 


fi2 ■ SCO, 3° 


c-2 

E ■ 9EBV 




E 2 


■ !EBV 11CE 


E 2 


' 9ia 


0.10-0.14 


1 


.35 


± 


07 


1 


09 


± 


11 


1 


20 


± 


11 


1 





±0.7 





1 ± 0.7 




9 ± 6 




7 ± 66 


0.00 ± 


03 





96 ± 0.18 


1.68 


± 11 


0.14-0.20 


1 


56 


± 


13 


1 


33 


± 


15 


1 


55 


± 


17 







±0.8 




2 ± 1.8 


7 


1 ± 3 


.0 


± 1 


0.79 ± 


33 





44 ± 0.29 


1.97 


± 


0.20-0.28 


1 


82 


± 


10 


1 


37 


± 11 


11 


1 


66 


± 0. 


13 


1 


2 


±0.6 


6 


2 ± 1.3 


6 


1 ± 2 





± 1 


1.11 ± 


23 





24 ± 0.20 


1.83 


± 11 


0.28-0.40 


2 


00 


± 0. 


09 


1 


70 


± 


09 


1 


69 


± 0. 


12 


1 


7 


±0.6 


9 


7 ± 1.2 




5 ± 1 


.6 


13 ± S 


1.19 ± 


19 





55 ± 0.17 


1.83 


± 


0.40-0.56 


1 


95 


± 


OS 


1 


76 


± 


09 


2 


11 


± (1 


11 


1 




±0.5 


8 


4 ± 1.1 


7 


± 1 


4 


9 ± 6 


1.33 ± 


1 1 





83 ±0.14 


1.68 


± 


0.56-0.80 


2 


10 


± 


OS 


1 


76 


± 


08 


2 


01 


± 


11 


2 


1 


±0.5 


7 


6 ± 1.0 


6 


4 ± 1 


3 


4 ± 5 


1.33 ± 


15 





92 ± 0.14 


1.42 


± 


0.80-1.13 


1 


90 


± 0. 


OS 


1 


70 


± 


08 


1 


95 


± 


10 


2 


4 


±0.5 


9 


6 ± 1.0 


6 


1 ± 1 


3 


17 ± 5 


1.03 ± 


14 





88 ± 0.14 


1.41 


± 


1.13-1.60 


1 


79 


± 0. 


0!) 


1 


59 


± 


09 


1 


90 


± 


10 


1 


1 


±0.5 


6 


7 ± 1.0 


5 


± 1 


3 


11 ± 5 


1.05 ± 


14 





79 ± 0.14 


1.23 


± 


1.60-2.26 


1 


74 


± 


09 


1 


36 


± 


09 


1 


15 


± 


11 


2 





±0.5 


6 


6 ± 1.0 


5 


± 1 


3 


7 ± 4 


0.98 ± 


15 





57 ± 0.14 


0.85 


± 


2.26-3.20 


1 


37 


± 


10 


1 


13 


± 


09 


1 


25 


± o 


11 


1 


5 


±0.5 


3 


7 ± 1.0 


3 


8 ± 1 


3 


3 ± 4 


0.80 ± 


11 





57 ± 0.14 


0.79 


± 11 


3.20-4.53 





84 


± 


10 





85 


± 


09 





99 


± 0. 


11 





1 


±0.5 


6 


1 ± 1.0 


4 


± 1 


3 


11 ± 4 


0.71 ± 


1 1 


0. 


43 ± 0.14 


0.91 


± 


4.53-6.40 





65 


± 


10 





48 


± 


09 





74 


± o 


11 





7 


±0.5 


3 


2 ± 1.0 


3 


9 ± 1 


3 


± 


0.65 ± 


15 


0. 


13 ± 0.15 


0.84 


± 


6.40-9.05 





54 


± 


10 





34 


± 


09 





37 


± 0. 


10 





9 


±0.5 


2 


5 ± 1.0 


1 


1 ± 1 


1 


2 ± 3 


0.55 ± 


14 





06 ± 0.15 


0.59 


± 


9.05-25.6 





38 


± (1 


OS 





23 


± 


13 





13 


± o 


10 





7 


± 0.4 


2 


3 ± 0.6 


1 


5 ± 


7 


3 ± 2 


0.1 ± 


3 




4 ± 1 


0.45 


± 



a Units: E 2 q HJ ^lO -24 McV 2 s' 1 sr -1 McV 1 ), E 2 - q co ^lO -4 McV 2 s _1 cm -2 sr -1 McV" 1 (K km s -1 ) -1 ), E 2 

g EBV (10 -2 MeV 2 s _1 cm -2 sr -1 McV -1 mag" 1 ), E 2 ■ q iao (10~ 3 McV 2 s _1 cm -2 sr -1 McV -1 ) 

The subscripts refer to four regions defined to perform the analysis; 1) Local arm, 2) interarm region, 3) Perseus arm, 4) beyond the Perseus arm. 
c Some of parameters are consistent with and thus are not well determined. Wo present them for reference. 
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Table 3 

Number of counts in each energy bin. 



Energy (GeV) Counts 






10- 


■0 


14 


26673 





14- 


■0 


20 


71637 





20- 





28 


91336 





28- 


■0 


40 


93286 





40- 





50 


78330 





56- 


■0 


80 


61337 





80- 


■1. 


13 


45386 


1. 


13- 


■1. 


00 


30713 


1. 


60- 


2 


20 


19351 


2 


26- 


3 


20 


11301 


3 


20- 


■4. 


53 


6426 


4 


53- 


6 


40 


3761 


6 


40- 


9 


05 


2095 





05- 


25.6 


2333 



Table 4 

Log-likelihood and emissivities for several choices of T s . 



Ts 


ln(L) 


9hi,i(E > 100 MeV) 


<Jhi,2(£ > 100 MeV) 


<7m,3(-E > 100 MeV) 


100 K 


114407.6 


1.32 ± 0.04 


1.27 ± 0.05 


0.86 ± 0.06 


125 K 


114480.1 


1.47 ± 0.05 


1.26 ± 0.06 


1.14 ± 0.08 


250 K 


114533.8 


1.62 ± 0.04 


1.35 ± 0.05 


1.53 ± 0.06 


400 K 


114544.5 


1.67 ± 0.07 


1.39 ±0.08 


1.64 ± 0.09 


optically-thin 


114552.8 


1.70 ± 0.07 


1.39 ±0.07 


1.77 ±0.09 



a Units: gHi,i(10 26 photons s 1 sr 1 H-atom 1 ) 

b Thc subscripts refer to four regions defined to perform the analysis; 1) Local arm, 2) interarm region, 3) Perseus arm. 
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